QUASIDIAGONALITY AND THE FINITE SECTION METHOD 
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Abstract. Quasidiagonal operators on a Hilbert space are a large and important class 
(containing all self-adjoint operators for instance). They are also perfectly suited for study 
via the finite section method (a particular Galerkin method). Indeed, the very definition of 
quasidiagonality yields finite sections with good convergence properties. Moreover, simple 
operator theory techniques yield estimates on certain rates of convergence. In the case of 
quasidiagonal band operators both the finite sections and rates of convergence are explicitly 
given. 



1. Introduction 

This paper follows up on an observation from [4, Section 10]. Namely we see what happens 
when one takes the general theorems from [9] and applies them to the particular class of 
quasidiagonal operators. The end results are quite satisfying both from a theoretical point 
of view and a practical computational point of view. We will see that the finite section 
method works very well on the class of quasidiagonal operators in general and some explicit 
and reasonable rates of convergence can be worked out in some interesting cases. 

Let H be a separable, complex Hilbert space, B(H) denote the bounded, linear operators 
on H and T 6 B(H) be given. An operator T is called quasidiagonal if there exist finite 
rank orthogonal projections Pi, P 2 , P3 . . . such that ||P„(t>) — v\\ ^ for all vectors v £ H 
and which asymptotically commute with T; 

\\[T,P n }\\ = \\TP n - P n T\\ -+ 0. 

This notion is due to Halmos (cf. [11], [3], [16]) and is inspired by the classical Weyl-von 
Neumann Theorem; if S G B(H) is a self-adjoint operator then S = D + k where D is a 
diagonal operator with respect to some orthonormal basis of H and A; is a compact oper- 
ator. Indeed, Halmos observed that T e B>(H) is quasidiagonal if and only if T = B + k 
where B is a block diagonal operator with respect to some orthonormal basis of H and 
k is compact. (This means there exist projections as in the definition of quasidiagonality 
but with the stronger properties that P n H C P n+ \H and [P n , B] — for all n - writing 
H = PxH © (P 2 - Pi)P © (P 3 - P 2 )H © • ■ ■ , the matrix of B with respect to this de- 
composition is easily seen to be block diagonal.) In particular, every self-adjoint operator is 
quasidiagonal and, more generally, a result of Berg implies that every normal operator is also 
quasidiagonal (cf. [8, Corollary II.4.2]). The self-adjoints contain a large set of interesting 
examples (Almost Mathieu operators, discretized Hamiltonians, etc.) and there are plenty 
of interesting non-self-adjoint, but still normal, examples like difference operators coming 
from PDE's with constant coefficients. There are even interesting non-normal, but still qua- 
sidiagonal, operators such as certain weighted shifts and Toeplitz operators (not all of these 
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are quasidiagonal but many examples are known to be quasidiagonal) to which the ideas of 
this paper apply. 

Section 2 of this note contains a rather extensive list of definitions and recalls some of the 
theorems from [9] that we will need. There is nothing new in this section but we hope the 
unfamiliar reader will find it a convenient reference. 

Sections 3, 4, 5 and 6 are devoted to some general theory. These sections are the theoretical 
foundation on which the actual numerical analysis will be based. That is, we make some 
very general observations which are of little computational value by themselves. In some 
sense, one can regard these results as analogues of an existence result for a large class of 
PDE's. When faced with an actual PDE it is, of course, good to know via some general 
theory that a solution exists and even better if one knows it is unique. However, the real 
goal of numerical analysis is to find explicit algorithms which are guaranteed to converge 
to a solution and, moreover, study the rate of convergence. With this analogy in mind one 
can regard sections 3 through 6 as proving existence and uniqueness results as well as giving 
general estimates on rates of convergence. Then in sections 7 and 8 some explicit algorithms 
are described which put the general machine into motion in two special cases. Below we give 
a more detailed description of the contents of this paper. 

In Section 3 we observe the existence of good finite sections for quasidiagonal operators 
and the uniqueness of the limits of things like norms, condition numbers, spectral quantities, 
etc. of the corresponding finite sections. 

Section 4 gives some simple estimates on rates of convergence for the finite section method. 
Though the rates of convergence are not always very good they are best possible (at this 
level of generality) as shown by simple examples. (However, in certain cases the rates of 
convergence can be quite good.) Section 5 simply takes the results of sections 3 and 4 
and specializes them to the cases of normal and self-adjoint operators where the rates of 
convergence have cleaner statements. In section 6 we quickly point out that the resolution 
of a well-known problem in operator theory (Herrero's approximation problem) yields an 
improved existence result for many quasidiagonal operators. Putting this result into practice 
in special cases seems quite difficult but a related strategy, inspired by the observations of 
section 6, will be the subject of another paper (cf. [5]). 

Sections 7 and 8 are devoted to explicit descriptions of how to implement the quasidi- 
agonal approach in the cases of unilateral and, respectively, bilateral band operators. In 
the unilateral case we use canonical finite sections and translate our previous results into 
computable error estimates. In the bilateral case a famous technique of Berg turns out to 
work beautifully and, again, computable rates of convergence are given. 

Acknowledgement. I thank Bill Arveson and Marc Rieffel for sparking my interest in numeri- 
cal analysis and, more generally, for providing more support and inspiration to me than they 
probably realize. I am also indebted to the authors of [9] and [2] as they have made a large 
body of work readily accessible to the C*-community. Finally, I thank my colleague Ludmil 
Zikatanov for numerous helpful conversations. 

2. Notation, definitions and known results 

This section is devoted to setting notation, recalling a number of definitions and stating 
some theorems from [9] that will be required. 

Let T e B(H) be a bounded operator on some complex Hilbert space H and Pi < P-i < ■ ■ ■ 
be a filtration (i.e. each P n = P* = P% is a finite rank, orthogonal projection, the range of 
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P n is contained in the range of P n+ \ and \\P n (v ) — v\\ — > for all vectors v G H). If the rank 
of P n is some integer k(n) then we will often - but not always - identify P n B(H)P n with the 
complex k(n) x k(n) matrices Mfc(„)(C). We will refer to the sequence of operators (P n TP n ) 
as the finite sections (of T with respect to {P n })- The finite section method is simply to 
try and study the (infinite dimensional) operator T via its finite sections with respect to a 
suitable filtration. For example, can one recover the spectrum of T as some sort of limit of 
the spectra of (the finite dimensional matrices) P n TP n l Of course, the choice of filtration 
is very important as it can easily happen that the finite sections do not approximate the 
original operator in any reasonable sense. A classic example is that of the unilateral shift S 
on l 2 (N). With respect to the canonical basis, the matrix of S is below. 

/0 • • \ 

1 ••• 

1 ••• 

1 ••• 

V : ■■■/ 

Cutting S to the canonical filtration yields nilpotent matrices and thus the spectra of the 
finite sections are all one point sets {0} which do not converge to the spectrum of S (which 
is the unit disc). 

Some operator theory definitions. 

Definition 2.1. Given an operator T G P>(H): 

(1) For e > 0, the e-pseudospectrum of T is 

a^{T) = {A G C : ||(A - T)~ 1 \\ > -}, 

e 

where A — T denotes the operator which maps v G B(H) to \v — Tv and we will adopt 
the convention that if an operator X G B(H) is not invertible then = oo. In 

particular, the usual spectrum a(T) is contained in a^ e \T) for every e > and, 
moreover, 

a(T) = f> (£) (T). 

£>0 

See http://web.comlab.ox.ac.uk/projects/pseudospectra for more, including 
software for computing pseudospectra. 

(2) The singular values of T are a 2 (T) = {A > : A 2 G ct(T*T)}. 

(3) We let C*(T) denote the unital, C*-algebra generated by T. It is the smallest unital 
C*-subalgebra of B(H) which contains T and can be realized as the norm closure of 
the set of all polynomials (where T° is defined to be the identity operator) in T and 
T*. A state on C*(T) is a linear functional <p : C*(T) -> C such that ||0|| = <p(I) = 1, 
where / is the identity operator on H . The set of all states is denoted S(C*(T)). 

(4) The spatial numerical range of T is 

SN(T) = {{Tv,v) : v G H, \\v\\ = 1} 
and the algebraic numerical range of T is 

AN(T) = {<j>(T) : G S(C*(T))}. 
Note that SN(T) C AN(T). 
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(5) T is called Moore-Penrose invertible if there exists an operator T + G B(H) such 
that TT+T = T, T + TT + = T+, (TT+)* = TT+ and (T+T)* = T+T. If such T+ 
exists then it is unique and is called the Moore-Penrose inverse of T. It turns out 
that if T+ exists then T+ G C*(T) (cf. [9, Theorem 2.15]) and hence, for an abstract 
C*-algebra the notion of Moore-Penrose invertible does not depend on the choice 
of faithful representation. Also, recall that T+ determines least square solutions to 
operator equations; if y G H then x = T + y is the least square solution of the equation 
Tx = y (i.e. x is the unique vector such that \\Tx — y\\ = mi{_\\Tu — y\\ '■ u G H} and 
||x|| = inf{||u|| : \\Tu — y\\ = mi{\\Tu — y\\ : u G -P}})- 

(6) If T is invertible then the condition number of T is 

cond{T) = IITIHIT- 1 !!, 

while the generalized condition number of a Moore-Penrose invertible operator T is 
the number ||T||||T+ 1| . 

We now need to recall some definitions motivated by the finite section method, but this 
will require introducing a bit more notation. If Pi < P 2 < • • ■ is a filtration of H and the 
rank of P n is k(n) then it is useful to build a C*-universe where finite sections can be treated 
as a single element. Namely we let 

IIM fc(n) (C) = {(x n ) : x n G M fc(n) (C),sup ||x n || < oo} 

and 

0M fc(n) (C) = {(x n ) G nM fc(n) (C) : lim ||x n || = 0} 

denote the l°° and, respectively, Co direct sum of the M/ c ( n )(C)'s (where we identify P n B(H)P n 
and Mfc(„)(C)). Note that ©M fc ( n )(C) is a closed, two-sided ideal in UM k ( n )(C) and hence 
there is a canonical quotient mapping ir : IIM fc ( n )(C) — > IIM fc ( n )(C)/ © M fc („)(C). This 
quotient algebra turns out to be the right place to study finite sections. Note that if T G 
B{H) then sup ||P n TP„|| < oo (actually, = \\T\\) and hence the sequence of finite sections 
(P n TP n ) can be regarded as a single element in nM^( n )(C). The real element of interest, 
however, turns out to be ir((P n TP n )). 

Definitions arising from the finite section method. 

Definition 2.2. Let T G B(H), Pi < P 2 < ■ ■ ■ be a filtration of H and (P n TP n ) be the 
corresponding finite sections. 

(1) (P n TP n ) is said to be applicable if there exists an integer n such that for every 
vector y G H and integer n > no there exists a unique solution x n to the equation 
P n TP n x n = y n and, moreover, the vectors x n converge (in norm) to a solution x of 
the equation Tx = y. 

(2) (P n TP n ) is said to be stable if there exists an integer n such that for every integer 
n > n the matrices P n TP n are invertible and sup n> „ ||(P n TP n ) _1 || < oo. A basic 
result of Polski states that (P n TP n ) is applicable if and only if T is invertible and 
(P n TP n ) is stable (cf. [9, Theorem 1.4]). 

(3) The finite sections (P n TP n ) are said to be fractal if for every subsequence n\ < n 2 < 
n 3 --- there exists a *-isomorphism a : C*(-K((P n TP nj ))) — >■ C*(7r((P n TP„))) such 
that a(n((P nj T B nj ))) = n((P n TP n )). (We have abused notation here in using n to 
denote both of the canonical quotient maps n n M fc ( n )(C) — > n n M fe („)(C)/ ©„M fe ( n )(C) 
and UjM k{nj) (C) -> n,M fe(n . ) (C)/ Af fc(n .)(C).) 
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(4) For a matrix S E M k (C) and e > 0, the e-regularization of S is denoted S £ and 
defined as follows. Let S = UDV be the singular value decomposition of S (i.e. 
U, V e Mfc(C) are unitaries and P = diag(Xi, . . . , A&) where Ai < A2 < • • • < X k 
are the (ordered) singular values of S. Define D £ = diag(Xi, . . . , X k ) where Xj = Xj 
if Xj > e and Xj = otherwise. The e-regularization of S is then defined to be 
S £ = UD £ V. 

(5) The finite sections (P n TP n ) are said to be stably regularizable if there exists a number 
e > such that || (P n TP n ) £ — P n TP n \\ — > and sup ||(P n TP n )+|| < 00. In this case 
a Polski-type theorem (cf. [9, Theorem 2.13]) states that for every vector y E H the 
least square solutions x n to the (regularized) equations (P n TP n ) £ x n = P n y converge 
in norm to the least square solutions of Tx = y (i.e. T is Moore-Penrose invertible 
and (P n TP n )fy — > T + y in norm for every y G H). 

One of the main topics of this paper is convergence aspects of spectra. As such it will be 
convenient to introduce some set notation and recall Hausdorff convergence of sets. 
Approximate inclusions and Hausdorff convergence of sets. 

Definition 2.3. Let C denote the complex plane. 

(1) For a set S C C and 5 > we will let Ns(H) denote the ^-neighborhood of E (i.e. 
the union of all 5-balls centered at points of E) and A^(E) be its closure. 

(2) Given two sets E, A C C we say E is 5-contained in A if E C Ns(A). We prefer (and 
will use) the following notation in this situation; E C 5 A. 

(3) Given two compact sets E, A C C their Hausdorff distance is 

eta (E, A) = maxjsup d(a, A), sup d(X, E)}, 
o-gs AeA 

where d(a, A) = infA G A \cr — X\. Note that 5 > du (E, A) if and only if E is 5-contained 
in A and vice versa - i.e. A C 5 E and E C 5 A. 

(4) Given a sequence of compact sets E n C C and another compact E C C we will write 

E n — > E 

if d H (E n , E) — > - i.e. if for each 5 > we have that A n C 5 A and A C 5 A n for all 
sufficiently large n. 

We close this section by stating all of the theorems from [9] that we will need. Roughly 
speaking these results say that all convergence questions for the finite section method have 
a C*-formulation in nM fc(n) (C)/ © M k{n) (C). 

Theorem 2.4. (Operator Equations) 

(1) (Kozak's Theorem) (P n TP n ) is stable if and only if n((P n TP n )) is an invertible e/e- 
ment in UM k{n) (C)/ © M fc(n) (C). 

(2) (P n TP n ) is stably regularizable if and only ifir((P n TP n )) is Moore-Penrose invertible 
in nil4 (n) (C)/©M fe(n) (C). 

See [9, Theorem 1.15] for a proof of the first statement above and [9, Theorem 2.19] for 
the second. 

Theorem 2.5. [9, Theorem 2.28] If (P n TP n ) is fractal and stably regularizable then the 
generalized condition number of n((P n TP n )) is equal to lim ||P n TP n || ||X n ||, where (X n ) is 
any sequence such that ir((X n )) is the Moore-Penrose inverse of ir((P n TP n )). 
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Theorem 2.6. (Spectral Approximations) If (P„TP n ) is fractal then the following hold: 

(1) For every e > 0, <r&(P n TP n ) -> ^ £ \tv( (P n TP n ))). 

(2) a 2 (P n TP n ) -> a 2 (7r((P n TP n )))- 

(3) If each P n TP n happens to be normal then a(P n TP n ) — > a(ii((P n TP n ))). 

The first statement above is a combination of Theorems 3.31, 3.33(b) and Proposition 3.6 
in [9]. The second is a consequence of Theorems 3.22, 3.23 and Proposition 3.6 in [9]. The 
last statement follows from Corollary 3.18, Theorem 3.20 and Proposition 3.6 in [9]. 

3. Existence and Uniqueness Results 

Let Pi < P2 < • • • be a filtration of H . If the rank of P n is some integer k(n) and we 
identify P n B(H)P n with the complex k(n) x k(n) matrices Mfc(„)(C) then we get a natural 
linear mapping <p n : B(H) — > M fc ( n )(C) defined by <p n {X) = P n XP n for all X G B(H). 
We further define linear maps $ : -> IIM fc(n) (C), $(X) = (</>„(X)) and p : B(H) ^ 

UM k{n) (C)/ © M fc(n) (C), p = 7T o $ where vr : IIM fc(n) (C) -> IIM fc(n) (C)/ © M fc(n) (C) is the 
canonical quotient mapping. 

The following lemma, which is well known to the operator algebra community, is really 
the key observation of this section. The proof is also very simple so we include it for the 
reader's convenience. 

Lemma 3.1. // || [T, P n ] || — >• then p restricts to an injective *-homomorphism on the unital 
C* '-algebra, denoted C*(T), generated byT. In other words, p\c*(T) '■ C*(T) — > L T Mfe( ri )(C)/ © 
Mk( n )(C) is a ^-isomorphism onto it's image. 

Proof. First note that the projections P n asymptotically commute with T* as well (since 
|| [T,P n ] || = || [T,P n ]* || = || [T*, P n ] ||). A simple interpolation argument together with the 
triangle inequality shows that (P n ) asymptotically commute with any polynomial in T and 
T*. It follows that the P n 's asymptotically commute with every element in C*(T). Since 
||p(X)|| = ||X|| for all X e B(H) (in particular, p is injective) it only remains to show that 
p restricts to a multiplicative map on C*(T). This follows easily since for all a, b G C*(T) 
we have \\(f> n (ab) - (f) n (a)(f) n (b)\\ = \\P n abP n - P n aP n bP n \\ = ||P n (P n a - aP n )6P„|| -> 0. □ 

Lemma 3.2. Let T e B(H) be a quasidiagonal operator and P± < Pi < • • • be any increasing 
sequence of finite rank projections such that ||[P„,T]|| — > and \\P n (v) — v\\ — > ; for all 
v G H , as n — > 00. Then (P„TP n ) is fractal. 

Proof Every subsequence of {P n } also asymptotically commutes with T and converges 
strongly to the identity. Hence the corresponding map p : C*{T) — > UM k ( n .)(C) / ©M fe („ 3 )(C) 
is a *-isomorphism. □ 

As stated in the title of this section, we regard the following results as existence and 
uniqueness results. That is, for a quasidiagonal operator one can always find a sequence of 
projections such that the finite sections converge where we want (existence) and, in fact, 
the finite sections arising from any other asymptotically commuting sequence will converge 
to the same thing (uniqueness). 

Theorem 3.3. (Operator Equations) Let T G B(H) be a quasidiagonal operator and Pi < 
P 2 < • • ■ be any filtration of H which asymptotically commutes with T. Then the following 
statements hold: 

(1) The sequence (P„TP n ) G IIMfc( n )(C) is stable if and only ifT is invertible. 
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(2) The sequence (P n TP n ) G IIM fe ( n )(C) is applicable if and only if T is invertible. 

(3) The sequence (P n TP n ) is stably regularizable if and only if T is Moore-Penrose in- 
vertible (i.e. the range ofT is closed - see [9, Theorem 2.4]/ 

(4) Assume that T has closed range, e > is some number strictly smaller than inf{A : 
< A G (T 2 (T)} 7 y G H is given and x G H is the least square solution to the equation 
Tx = y. Ify n = P n y and x n is the least square solution to (the e-regularized equation) 
(P n TP n ) £ x n = y n then \\x n — x\\ — > as n — > oo. 

Proof. (1) follows from Kozak's Theorem and Lemma 3.1 since p(T) is invertible if and only 
if T is invertible. 

(2) follows immediately from (1) and Polski's Theorem which was stated right after the 
definition of stability in the previous section. 

(3) follows from the second part of Theorem 2.4 and Lemma 3.1. 

(4) follows from (3) and the least square version of Polski's Theorem stated after the 
definition of stably regularizable. Strictly speaking, [9, Theorem 2.13] only ensures the 
existence of an e > for which our claim holds true. However, the convergence of spectra 
results below will ensure that any < e < inf{A : < A G a(T*T)} will do. □ 

Theorem 3.4. (Norms and Condition Numbers) Let T G B(H) be a quasidiagonal operator 
and Pi < P 2 < • • ■ be any filtration of H which asymptotically commutes with T . Then the 
following statements hold: 

(1) ||T|| = lim \\P n TP n \\ and \\T~ l \\ = lim ||(P„TP n ) _1 ||- 

(2) IfT is invertible then cond(T) = lim cond(P n TP n ). 

71-+00 

(3) IfT has closed range, e > is some number strictly smaller than inf{A : < A G 
o"2(T)} ; and T + (resp. (P n TP n )f) denotes the Moore-Penrose inverse of T (resp. 
(P n TP n ) £ ) then the generalized condition numbers of (P n TP n ) £ converge to the gen- 
eralized condition numbers ofT - i.e. ||T||||T + || = lim ||(P n TP„) £ || ||(P„TP„)+||. 

Proof. The first part of statement (1) is trivial. In fact, for every X G B(H), the norms of 
P n XP n increase up to \\X\\ (since P n P n+ i = P n and thus ||P n XP„|| = \\P n (P n+ iXP n+1 )P n \\ < 
||P n+ iXP n+ i||). For the inverse statement we first assume that T is invertible and note that 
IIP- 1 ! = lim||P n T- 1 P„||. Also, since T _1 G C*(T), the fact that ptT" 1 ) = p(T)" 1 implies 
that {PnT^P^-^T^P^' 1 ) G ©M fc(n) (C). It follows that \\PnT~ 1 P n - {P n T~ l P^W -> 
as n -> oo and hence \\T- l \\ = lim \\P n T^ l P n \\ =]im\\(P n TP n )- x \\. If T is not invertible but 
lim ||(P n TP n ) _1 || 7^ oo then we can find a subsequence {rij} such that lim ||(P n TP„ 7 ) _1 || = 

n^oo j— *oo 

c < oo. However, the projections {P U] } also asymptotically commute with T and the 
finite sections (P nj TP n .) are now stable. Thus, by Theorem 3.3 (1), T is invertible - a 
contradiction. 

(2) follows immediately from (1). 

(3) follows from Theorem 2.5 and Lemmas 3.1 and 3.2. □ 

Theorem 3.5. (Spectral Approximations) Let T G B(H) be a quasidiagonal operator and 
Pi < Pi < • • • be any filtration of H which asymptotically commutes with T . Then the 
following statements hold: 

(1) For every e > 0, a^(P n TP n ) -> a^(T). 

(2) a(T) = H (lima( £ )(P n TP„)^, where lim (P n TP„) = a^(T) by (1). 
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(3) a 2 (P n TP n )^a 2 (T). 

Proof. Given Lemmas 3.1 and 3.2, the proof is reduced to recalling Theorem 2.6 from the 
last section. □ 

Remark 3.6. As mentioned earlier, a theorem of Berg states that every normal operator on 
a Hilbert space is quasidiagonal and hence the existence and uniqueness results above can 
be applied. Even in the case of a self-adjoint operator the existence of a filtration for which 
the finite sections are guaranteed to converge nicely appears to be a new result. 

Remark 3.7. Many standard examples of operators are not quasidiagonal. For example, 
any Fredholm operator with non-zero index is not quasidiagonal. However, there are very 
nice characterizations of quasidiagonality for essentially normal operators (cf. [8, Corollary 
IX.7.4, Section LX.8]) and weighted shift operators ([14], [13]). 

Remark 3.8. Though many operators are not quasidiagonal, a remarkable result of Voiculescu 
implies that if one starts with an arbitrary T G B(H) then it is possible to construct 
a quasidiagonal operator out of T. Namely, let Ai, A 2 , . . . be any listing of all the rational 
numbers in the half open interval (0, 1]. If T G B(H) is arbitrary then Voiculescu's homotopy 
invariance theorem (cf. [16]) implies that 

T = (BienKT, 
acting on (B ieN H, is a quasidiagonal operator. 

We close this section with a result which has nothing to do with quasidiagonality. Namely, 
we will consider the convergence of the algebraic numerical ranges of finite sections. It is well 
known that spatial numerical ranges are well behaved (cf. [9, Theorem 3.52]), but algebraic 
numerical ranges also behave nicely. Though we would be a bit surprised if this result is not 
already known, it is not mentioned in [9] and hence we include a proof. 

Theorem 3.9. (Numerical Ranges) Let T G B(H) be arbitrary and Pi < P 2 < . . . be any 

sequence of finite rank projections such that \\P n {v) — v\\ — > 0, as n — > oo ; for all v G H. 
Then the following statements hold: 

(1) For every n, SN(P n TP n ) = AN(P n TP n ) C SN(T ) C AN(T). 

(2) SN(P n TP n ) = AN(P n TP n ) -> AN(T) = SN(T) = AN(T). 

Proof. It is immediate from the definitions that SN(X) C AN(X) for any X G B{H). 
Moreover, it is also known that AN(X) = SN(X) (cf. [9, Theorem 3.45]). Since the unit 
ball of a finite dimensional Hilbert space (e.g. P n H) is (norm) compact, these two facts imply 
that SN(P n TP n ) = AN(P n TP n ) (for a given state 0, choose unit vectors Vk G P n H such 
that (j)(P n TP n ) = lim(P n TP n Vk,Vk) and any cluster point of {vk} will do the trick). The 
proof of (1) is then complete after we observe that SN(P n TP n ) C SN(T): for a unit vector 
v G P n H we have (Tv,v) = (TP n v,P n v) = (P n TP n v,v). 

Statement (2) follows from (1) and the general fact that SN{P n TP n ) -> SN{T) (cf. [9, 
Theorem 3.52]). □ 

4. Commutators and Rates of Convergence 

In this section we observe that norms of the commutators ||[P n ,T]|| provide a simple and 
natural upper bound for the rate of convergence of finite sections (P n TP n ). For spectral 
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approximations this only yields 'one sided' rates of convergence but, sadly a trivial example 
at the end of this section shows that this is best possible. 

From a computational point of view this will, in general, not provide a very useful measure- 
ment. However, we will see in later sections that in some cases these commutator estimates 
yield computationally relevant rates of convergence. 

As was the case for Lemma 3.1, the following (well known) lemma is really the key remark 
of this section. 

Lemma 4.1. Let T G B(H) be arbitrary and P G B(H) be a self-adjoint projection. Then 

\\T -(PTP + P L TP L )\\ < \\PT-TP\\, 

where P 1 - = I — P. 

Proof. Since T = PTP + PTP ± + P ± TP + P ± TP ± (T = IT I and / = P + P^) we have that 
T — (PTP + P^TP^) = PTP 1 - + P^TP. Since the domains and ranges of the operators 
PTP ± and P L TP are orthogonal we have \\PTP ± + P L TP\\ = max.{\\PTP ± \\, || P^TPH}. 
However, each of the latter norms is dominated by || PT — TP\\ - for example, HP^TPH = 
\\P ± TP-P ± PT\\<\\P ± \\\\TP-PT\\. □ 

Though we prefer to state the results below in terms of norms of commutators, we should 
also point out that these norms are determined by the norms of finite dimensional matrices: 

\\PT-TP\\ = \\PTP ± - P^TPW 

= max{||PTP ± ||, ||P- L TP||} 

= max{||P- L r*P||,||P- L TP||} 

= max{ || PTP^T'P || 1/2 , WP^P^TPf 12 }. 

Proposition 4.2. (Operator Equations) Let T G B(H) be a quasidiagonal operator and 
Pi < P2 < • • • be any filtration such that \\[P n ,T]\\ — > 0. Then the following statements hold: 

(1) IfT is invertible and ||[P n ,T]|| < p^rij then P n TP n is also invertible. 

(2) IfT is invertible, \\[P n ,T]\\ < y eH is given, x G H is the unique solution to 
Tx = y, y n = P n y and x n is the unique solution to P n TP n x n = y n then 

\\x-x n \\ < WT-'lKWy - y n \\ + ||x n ||||[P n ,T]||). 

Proof. The first statement follows from the following inequalities: 

\\P n TP n T- 1 P n -P n \\ = \\P n (TP n -P n T)T- 1 P n \\ < HtP^TlllIlT- 1 !! < 1. 

Indeed, this shows that P n TP n T- x P n = (PnTP^iPnT^Pn) is invertible (we now regard P n 
as the unit of the matrix algebra P n B(H)P n ) and hence P n TP n is right invertible. However, 
in finite dimensions right invertible implies invertible. 
The second assertion is also straightforward: 

||x n -x|| = \T~ x (Tx n -Tx)\ 

< ||T _1 ||||Tx n - y\ 

|| P || ||PP n x n P n TP n x n -\- y n y\\ 

< \\T^\\ (\\y n - 2/|| + ||x n || ||(TP n - P„T)P n ||) 

< ||T- 1 ||(||y-i frl || + || a ; B ||||[P B ,7l||). 

□ 
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Remark 4.3. (Computability) We regard the estimates above as reasonably computable ob- 
jects (at least in concrete cases). Indeed, if T and {P n } are concretely given then ||[P n ,T]|| 
is computable in terms of two finite dimensional matrices. (Warning: In this remark we 
take the point of view that anything in finite dimensions is 'computable' and limits of fi- 
nite dimensional calculations are regarded as 'reasonably computable'.) Also, note that 
\\y — P n y\\ = a/IMI 2 — ||-Pny|| 2 is computable (at least if one knows ||y||). 
Recall that the polar decomposition implies that for any invertible operator S we have 

||S _1 || = sup I 

Hence the quantities 1 1 CZ - ' 1 1 1 in the estimates above are reasonably computable since we 
know the singular values of P n TP n converge to the singular values of T. Also, < 
||(P n TP n ) _1 || \\y n \\ and thus can be estimated in concrete cases. 

Remark 4.4. (Least Square Estimates) We have been unable to find any 'reasonable' estimates 
on the rate of convergence of least square solutions. Our feeling is that this may be possible, 
however we don't yet see how to avoid the Moore-Penrose projection (a certain spectral 
projection for T*T - see [9, Theorem 2.15]) which arises in the construction of T + . 

Proposition 4.5. (Spectral Approximations) LetT G B(H) be a quasidiagonal operator and 
Pi < P2 < • • • be any filtration such that \\ [P n , T]\\ — > 0. Then the following statements hold: 

(1) For every e > 0, a^(P n TP n ) C g^+\ \[Pn,T\\\)(T). 

(2) a 2 (P n TP n ) C s - a 2 (T), where 5 n = y/2\\T\\ || [P„, T] || . 

Proof. We will need the following characterization of e-pseudospectra (cf. [9, Theorem 3.27]): 
If X G P(P) then 

^\X)= |J a(X + Y). 

YeB(H), \\Y\\<e 

Given this fact, it is clear that if \\X-Z\\ < 5 then a^(X) C a^^Z). Since a^(P n TP n ) C 
o~( e \P n TP n + P^TPn), statement (1) follows from Lemma 4.1 and the remarks above. 

Part (2) requires the following well known fact: If X, Y G B(H), X is self-adjoint (even 
normal will suffice) and ||-X"— Y \\ < S then <r(Y) C s o~(X). To complete the proof we first note 
that a(P n T*P n TP n ) C a(P n T*P n TP n + P^T*P^TP^) = a((P n T*P n + P^T*P^)(P n TP n + 
P^TP^)). However, Lemma 4.1 implies that ||(P n T*P n + P^T*P^)(P n TP n + P^TP^) - 
T*T\\ < 2||T|| || [P n ,T] || . Combining all these facts we see that a(P n T*P n TP n ) c 2 l™[ p "' T ]" 
a(T*T) - i.e. if A G (72 (P„TP n ) then there exists some point 7 G 02 (P) such that |A 2 — 7 2 | < 
2||T|| || [P n , T] \\. Since (A — 7) 2 < |A 2 — 7 2 |, statement (2) follows. Note that for large singular 
values one does not have to take square roots, however the estimate (A — 7) 2 < |A 2 — 7 2 | 
becomes sharp at zero. □ 

We end this section with a simple example which shows that there is no hope of giving 
rates of convergence of the quantities not mentioned in this section in terms of norms of 
commutators. 

Let j\ < j'2 < J3 • • • be integers and consider the diagonal, self-adjoint operator on Z 2 (N) 
whose matrix (w.r.t. the canonical orthonormal basis) is 



diag(2, . . . , 2, 1 + 1/2, 3 - 1/2, 2, . . . , 2, 1 + 1/3, 3 - 1/3, 2, . . . , 2, 1 + 1/4, 3 - 1/4, 2 • • • ), 
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where the first string of 2's has length j ± , the second string of 2's has length j 2 and so on. 
One should have in mind here that the jVs are growing like exponentials (or faster, if you 
like) in k. Now let P n denote the orthogonal projection onto the span of the first n basis 
vectors of l 2 (N). 

Note that [P n ,T] = for all n. On the other hand, the norms of P n TP n and it's inverse 
approach the norms of T and it's inverse as slowly (in n) as you like. Also, the condition 
numbers of P n TP n are approaching the condition number of T as slowly as you like. 

It would also be nice if we could approximately reverse the inclusions in Proposition 4.5 - 
i.e. if we knew something like o~(T) C a^ e \P n TP n ) for all n larger than some number which 
depended on commutators. This would amount to saying how quickly the pseudo-spectra 
of P n TP n "fill up" the pseudo-spectrum of T. Unfortunately, the example above shows that 
this can happen as slowly as desired. Similarly the numerical ranges of P n TP n will "fill up" 
the numerical range of T at as slow a rate as one wishes to prescribe. 

Finally, it would be nice for solving operator equations if we could somehow relate \\y — y n \\ 
and || [P n ,T] || but it is not hard to see that this is also impossible to do in general. 

5. Normal and Self-Adjoint Operators 

Here we observe that some of the statements from the previous section can be cleaned up 
a bit in the case that T is a normal or, better yet, self-adjoint operator. At the end of this 
section we also state a result from [4] which provides a very general existence result which is 
analogous to the Szego-type results obtained by Bill Arveson for certain self-adjoints (cf. [9, 
Chapter 7]). 

Proposition 5.1. If T e B(H) is a normal operator and P\ < Pi < ■ • • is any filtration 
such that \\[P n ,T]\\ — > then: 

(1) For every e > 0, cr {e \P n TP n ) -> N e (a(T)). 

(2) For every e > 0, <T&(P n TP n ) C Sn a(T), where S n = e + ||[T,P n ]||. 

Proof. For a normal operator T we have N e (o~(T)) = a^ e \T). This follows from the fact that 
if T is normal and A G C then || (A — T) _1 1| = d( - A 1{t)) ' wnere d(\, &(T)) is the distance from 
A to er(T). Hence we have just reformulated the results of sections 2 and 3. □ 

Remark 5.2. On page 116 in [9] a general conjecture concerning convergence of spectra of 
'essentially normal' approximation sequences is formulated. Though we have been unable to 
resolve this conjecture we do wish to point out that an affirmative answer would lead to a 
nice improvement in the previous result. Namely, it would follow that a{P n TP n ) — > o~(T) - 
i.e. we wouldn't have to bother with pseudospectra in the result above. 

Proposition 5.3. IfT e B(H) is a self-adjoint operator and P± < P 2 < ■ • • is any filtration 
such that || [P ni T] \\ —> then: 

(1) a(P n TP n )^a(T). 

(2) a(P n TP n ) C 5 " a(T), where 5 n = \\[T,P n ]\\. 

Proof. Assertion (1) follows from part (3) of Theorem 2.6 and Lemmas 3.1 and 3.2. The 
second part follows from Lemma 4.1, the fact that a(P n TP n ) C a(P n TP n + P^TP^) and 
the first sentence in the proof of part (2) of Proposition 4.5. □ 
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Remark 5.4. Following Arveson's ideas, it is possible to also compute the essential spectrum 
of a self-adjoint operator via quasidiagonal finite sections (cf. [9, Chapter 7]). It does not 
seem easy, however, to develop any rate of convergence estimates in this case. 

We conclude this section with an existence result which is analogous to Szego's celebrated 
limit theorem concerning spectral distributions of self- adjoint Toeplitz operators. In [4] we 
studied various notions of invariant means on C*-algebras and observed in Theorem 10.4 
that 'quasidiagonal' invariant means can be recovered following the ideas first laid out by 
Bill Arveson (as explained in [9, Section 7.2]). To avoid technicalities we will just restate 
the result in the case of normal operators. 

Theorem 5.5. Let T e B(H) be a normal operator. Then there exists a filtration Pi < P 2 < 
• • • such that \\[P n ,T}\\ —> and with the property that for every Borel, probability measure 
H supported on the essential spectrum of T there exists a subsequence {P nk } such that 

tr(P nk ■ P nk ) -> n 

in the weak—* topology, where tr(-) is the (unique) tracial state on P nk B(H)P nk and tr(P nk ■ 
P nk ) is the state on C*(T) defined by x i— > tr(P nk xP nk ) . 

Proof. Since normal operators generate exact C*-algebras (even nuclear ones) and every trace 
on an abelian C*-algebra is a quasidiagonal invariant mean (cf. [4, Lemma 6.4]) this result 
is a special case of [4, Theorem 10.4]. □ 

6. Improved Existence Result when C*(T) is Exact 

Mainly due to the pioneering work of Kirchberg, exact C*-algebras have been at the fore- 
front of operator algebra theory for the past decade. (See [12], [18].) This is a technical 
notion which we do not wish to discuss here, however it turns out that most (but not all) of 
the basic examples of C*-algebras enjoy this property (commutative algebras, AF algebras, 
irrational rotation algebras, Cuntz algebras and all subalgebras of these examples). Similarly, 
it turns out that many standard examples of operators on Hilbert space generate C*-algebras 
with this property. For example, C*(T) is exact whenever T e B{H) is essentially normal 
(e.g. self-adjoint, unitary or normal) or if T is a band operator (either unilateral or bilateral). 
Moreover, it follows from Arveson's extension theorem that this class of operators is closed 
under norm limits and hence band dominated operators (i.e. norm limits of band operators) 
also generate exact C*-algebras. 

In the early 80 's Herrero asked whether every quasidiagonal operator T was the norm limit 
of operators T n such that C*(T n ) was finite dimensional for each n. This turns out to be 
false (cf. [15], [17]) but the precise description of operators with this approximation property 
was discovered a couple years ago (cf. [6], [7]). 

Theorem 6.1. Let T e B(H) be given. There exist operators T n e B(H) such that \\T — 
T n \\ — > and C*(T n ) is finite dimensional (for all n) if and only if T is quasidiagonal and 
C*(T) is exact. 

This result is relevant to this paper because of the simple structure of finite dimensional 
C*-algebras: they are just finite direct sums of finite dimensional matrix algebras. So if T 
is close to T n and C*(T n ) is finite dimensional then one immediately gets a matrix whose 
£-pseudospectrum, for example, approximates that of T (since we may regard T n as a direct 
sum of finite dimensional matrices when viewed as an element of C*(T n )). In this setting 
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rates of convergence are controlled by \\T n — T\\ and estimates on the Hausdorff distance 
of spectra (in the normal case, at least) can be given. Moreover, standard representation 
theory of finite dimensional C*-algebras gives a decomposition of the underlying Hilbert 
space as a direct sum of finite dimensional reducing subspaces for C*(T n ) and thus a way 
of getting appropriate finite sections. Of course, actually finding these approximating T n 's 
and corresponding finite sections in specific examples does not appear to be easy and hence 
we have little more to say about it at this point. However, an effective means of computing 
spectra, inspired by the observations of this section, in irrational rotation algebras is worked 
out in [5]. 

7. QUASIDIAGONAL UNILATERAL BAND OPERATORS 

Though the results of the previous sections are quite satisfactory from a theoretical point 
of view some practical numerical analytic questions have went unaddressed. For example, 
if an explicitly given operator is known to be quasidiagonal how does one explicitly write 
down asymptotically commuting finite rank projections? Can one control the ranks of these 
projections? In the generality of the previous sections these questions cannot be answered 
(even for a general self-adjoint operator this is probably hopeless). However, in special cases 
some progress can be made and this is the subject of the final two sections of this paper. 

In this section we consider certain band operators. As mentioned before, not all band 
operators are quasidiagonal (e.g. unilateral shift) but a fairly broad class of band operators, 
described below, can be attacked with this technology. 

Our Hilbert space will now be / 2 (N) with canonical orthonormal basis {5j}« e N- We will 
regard /°°(N) as acting on l 2 (N) by multiplication operators (i.e. / G Z°°(N) acts on l 2 (N) by 
sending <5; to f(i)5i). T G B(l 2 (N)) will be a weighted shift operator with weight sequence 
CK2, . . .) - i.e. (ai, a 2 , • • •) is a bounded sequence of complex numbers and T(Si) = aA+i- 
Our basic assumption is: 

Assume the existence of a subsequence {ik}km suc h that \a ik \ — > as k — > oo. 

Note that a unilateral weighted shift is quasidiagonal if and only if the assumption above 
holds. Below we will observe that this assumption implies quasidiagonality while the absence 
of such a subsequence implies that the weighted shift is a Fredholm operator with index — 1 
and, hence, is not quasidiagonal. 

The band operators considered will be of the form 

S = f- m (T*r + /- m +i(TT^ + • • • f-iT* + f + fiT + • • • + / m _!T— 1 + f m T m , 

where fj G Z°°(N), — m < j < m, and m is a positive integer. Note that S is a band operator 
whose main diagonal is just /q. The off diagonal entries of S are somewhat complicated 
to describe, but the 'coefficients' fj are arbitrary and hence operators of the form above 
are still are rather large class. The reader may find it useful to consider a weight sequence 
of only zeroes and ones to get a feel for the kinds of operators arising this way. In fact, 
any band operator whose matrix has "butterflies" tending to zero can be handled with the 
quasidiagonal approach. In other words, any band operator which is a compact perturbation 
of a block diagonal band operator can be handled this way but rates of convergence can be 
complicated to state precisely so we stick to the situation described above. 

For ease of notation we will let q(x) = x m + x m ~ l + • • • + x and q' be the derivative of 
this polynomial. We will also let M = rnax_ m <j< m {||/j||}. Finally, we define finite rank 
projections Pt as the orthogonal projections onto the span of {5±, . . . , 5i k }. 
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Theorem 7.1. Let T, S and {Pk} be as above. Then the following statements hold: 

(1) If S is invertible and \a ik \ < ^s^JM^rfi ^ en P^SPk is invertible. 

(2) If S is invertible, \ai k \ < pprp^p^ , V £ / 2 (N) is given, x G / 2 (N) is the unique 
solution to the equation Sx = y, Uk = PkU and Xk is the unique solution to PkSPkXk = 
y k then 

\\x k - x\\ < US- 1 )! (\\y - Vn\\ + ||z„||M</(||P||)K|). 

(3) For every e > 0, the e -pseudo spectra of P k SPk converge to the e -pseudo spectrum of 
S andaV(P n SP n ) C(T (s+W(l|T||)K 

(4) The singular values of PkSPk converge to the singular values of S and <7 2 (P n TP n ) C Sn 
a 2 (T), where 5 n = y/2\\T\\Mq'(\\T\\)\a ik \. 

Proof. First note that ||[Pfc,T]|| = ||[Pfc,T*]|| = \cti k \ — > as k — > oo. Using the estimates 

l|[P^]|| < iipr^-rpr^ii + iirpr^- 1 -^^!! < ||[,p, r]|| ||r|| ||[,p, r-?- 1 !!, a simple 

induction argument shows that 

\\[P k ,T j ]\\ < jWTnp^TjW = jlirii^iaij. 

Since S is a band operator whose lower diagonals are determined by the powers of T and 
upper diagonals are determined by the powers of T*, we may use the fact that ||[P, S]\\ = 
max{ || PSP 1 - 1|, II P ± SP ||} to arrive at the following inequality: 

||[P fc ,S]||<Mg'(||T||)K|. 

With this inequality in hand, the various statements above follow from the corresponding 
results in sections 3 and 4. □ 

Remark 7.2. Note that the rank of Pk is %k and hence PkSPk is an %k x %k matrix. Thus the 
numerical efficiency in this setting is directly tied to the rate of growth of {ik} and the rate 
of decay of {a ik } (as functions of k). 

Remark 7.3. Evidently the theorem above can be extended with no difficulty to the case 
of bilateral shifts on Z 2 (Z) where one assumes the existence of subsequences tending to zero 
both in the positive and negative directions. 

Remark 7.4. Of course, the numerical ranges of PkSPk also converge to the numerical range 
of S since the proof above shows that ||[Pfc, S]\\ — > 0. Note also that a similar result can 
be worked out for any element in C*(T) and (complicated) error estimates can be explicitly 
given for any polynomial in T and T*. 

Remark 7.5. Let / G C(T) be a continuous function on the unit circle and Tf be the 
corresponding Toeplitz operator on / 2 (N). Of course, Tf is rarely a band operator and 
the diagonals of the matrix of Tf never have subsequences tending to zero (since they are 
constant). However, it follows from BDF theory that Tf is quasidiagonal if and only if the 
winding number of the curve /(T) is zero about every point in the complement of /(T) (cf. 
[8, Sections V.l and IX. 7]). In this case there are asymptotically commuting finite rank 
projections out there and we feel it would be interesting (but perhaps hard?) to find them. 

We conclude this section with an illustration of the usefulness of this method. Namely, 
we will give a new proof of the fact that the spectra of quasidiagonal weighted shifts is as 
large as possible. It is well known that if T is any weighted shift and p(T) is the spectral 
radius of T then a(T) = {A G C : |A| < p(T)}. While our technology will only handle 
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the quasidiagonal case, the argument is simple and cute so we include a proof. It will be 
convenient to note that 



P(T) 



lim \\T l \\ 1/l 



lim 



sup \a m a m+1 

men 



l/l 



whenever T is a weighted shift with weight sequence (a±, ct2, • • •)• 

Proposition 7.6. Let T be the weighted shift described at the beginning of this section. Then 
a(T) = {A G C : |A| < p(T)}. 

Proof. Evidently we may assume that ||T|| < 1 by rescaling, if necessary. 

Since the spectrum of T can be no larger than claimed, it suffices, thanks to part (2) of 
Theorem 3.5, to show that if 



|A| < lim 

1^00 



sup \a m a m+ i ■ ■■ctm+i 

mGN 



then for every e > 0, A belongs to the e-pseudospectra of PkTPk for all sufficiently large k. 
In fact, we will show a bit more (as it requires no extra effort and simplifies notation). Let 
P n be the projection onto the span of the first n basis vectors of / 2 (N) and we will show that 
A belongs to the e-pseudospectra of P n TP n for all sufficiently large n. 

First note that for all non-zero A the matrix A — P n TP n is invertible. A straightforward 
computation shows that its inverse is given by the following matrix. 



1 i 

A 2 
aiot2 
A3 

A4 



nr 




1 

A 

0.1 

A 2 
A3 



nr 1 - 







1 

A 

03 

A 2 



nn—1 
3 • 








1 

A 



IL a 

A™" 3 









"11-1 
A 2 



o\ 







Our goal is to show that if < |A| < lim/^oo I sup meN \a m a m+ i 



a 



m+l 



l/l 



< 1 then for 



every e > the norms of the matrices above are eventually bigger than l/e. To prove this 
it suffices to show that the modulus of some entry in the matrices above is greater than l/e. 
This, however, is the case because if |A| < \a m a m+ i ■ ■ ■ a m+ i — 5 then 

/-, , r w \ot m ot m +i ' ' ' ct m +i\ \a m a m+ i ■ ■ ■ a m+ i\ 
(l+o) < —, < 



IAI 



|AI Z+1 



Hence, if e > is given, we choose / big enough that (1 + 5) 1 > l/e and then we will be able 
to find a string a m , . . . , a m+ i such that 



1 \a m a m+ i ■ ■ ■ a m+ i \ 
e 



\i+i 



However, the numbers 



W 1 



are entries in all of the matrices above (for sufficiently large n). □ 

Remark 7.7. Note that the proof above adapts easily to bi-lateral weighted shifts which have 
subsequences tending to zero in both the positive and negative directions. 
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Remark 7.8. If T denotes the bilateral shift on Z 2 (Z) and P n denotes the orthogonal projection 
onto the span of {5_ n , . . . , 5 n } C Z 2 (Z) then the finite sections P n TP n are all nilpotent 
matrices and hence have spectrum {0} for all n (which, in particular, do not converge to the 
spectrum of T). Moreover, the proof above shows that for every e > 0, liminf a^(P n TP n ) 
contains the entire unit disc (see [9, Definition 3.1]). Hence we can't recover the spectrum 
of T (i.e. the unit circle) by intersecting over e > as in part (2) of Theorem 3.5. 

8. Berg's Technique and Bilateral Band Operators 

In our final section we consider bilateral band operators and apply a technique of Berg to 
produce explicit finite sections with the right quasidiagonality requirements. Berg's technique 
is a well-known, and highly successful, way of taking a shift operator and approximating it 
by cyclic shifts of a shorter length. We will not attempt to explain the philosophy behind 
this technique (as this is well done in [8, Section VI. 4], for example) but rather will just 
write down how it goes and check that it works. The diligent reader who actually checks all 
of the details will likely see what is going on in the process. 

In [14] it was shown that a bilateral weighted shift is quasidiagonal if and only if the weight 
sequence is "block balanced". The definition below is nothing but the obvious extension of 
Smucker's notion to subsets of Z°°(Z). 

Definition 8.1. A set S C Z°°(Z) is called block balanced if for each e > and k G N there 
exist integers m < and n > such that m + k < and for each sequence (aj) G S we have 
\a m+j - a n+j \ < e for < j < k. 

This definition just says that we can find strings of arbitrary length which are arbitrarily 
close to each other in the positive and negative directions. 

We now describe the types of bilateral band operators to which the quasidiagonal approach 
will apply. Let U be the bilateral shift on / 2 (Z) (i.e. U(Si) = 5i + i for the canonical basis 
{o~i}iez of Z 2 (Z)). Assume {f- p , ■ ■ ■ , f-i, fo, fi, ■ ■ ■ , f P } C /°°(Z) is a block balanced set and 
consider the band operator 

i=- P 

Since we have assumed {f- p , ■ ■ ■ , f-i, fo, fi, ■ ■ ■ , fp} to be block balanced we can, for any 
sequence Sk — > and natural numbers Sk — ■> oo, find integers mk,rik such that + Sk < 
< n k and 

\fi(m k +j) - fi(n k +j)\ < e k , 

for all — p < I < p and < j < s k . Define constants a^ k = |(1 — exp ^-jpf ) and /3j jk = 
|(1 + exp ) f° r — J — s k- A calculation shows that the set of vectors {a ^5 mk + 
Po,k^n k , • • • , a Skjk 5 mk+Sk + P Skjk 5 nk+Sk } is an orthonormal set. Hence the set 

V k = {Si : m k + s k < % < n k } U {a 0:k 5 mk + (3o yk 5 nk) . . . , a Skyk 5 mk+Sk + f3 Sk:k 8 nk+Sk } 

is also orthonormal. 

Though it is not obvious, it turns out that the orthogonal projection onto the subspace 
spanned by V k almost commutes with T. We will denote this projection by P k . 
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Lemma 8.2. Keeping notation as above we have 



\\[T,Pk]\\< 



(2p+l)e k , 7T 



+ 



s k + 2 



i=- P 



Proof. Our first task will be to write down the 5x5 block matrix of P k with respect to 
the decomposition Z 2 (Z) = H\ © H 2 ® H 3 ® H 4 ® H 5 , where Hi = span{5j : j < rrik}, 
H 2 = span{5j : m k < j < m k + s k }, H 3 = span{5j : m k + s k < j < n k }, H 4 = span{5j : 
n k < j <n k + s k } and H 5 = span{5j :n k + s k < j}. 

A tedious but straightforward calculation shows that in this decomposition P k becomes 



/0 






where / is the identity matrix, 



' o' l( 

A k = 





A k 


ci 













/ 








c k 



B k 




V o 







B k 



and 



C k — 







V 












o\ 





V 



/IA)I 





V o 






101 








where f3j is the complex conjugate of f3j. Since the matrix of each /; is also block diagonal 
with respect to the above decomposition we may write 



fi 



(F[ 








F\ 


















Fi 



0\ 




V 



Fi 



Hence 
\\[Pk,j, 



maxH 



\F 2 C k 



F 4 C k \ 



\FlCl 



F\Ct\ 



< e fc max{|/5 a |, • • • , lA^J} < y- 



Unfortunately, the matrix of the bilateral shift U is a bit messy with respect to the 
decomposition we have used above. However, conjugating by U (i.e. considering U ■ U*) just 
has the affect of shifting the matrix down the diagonal one entry. This remark together with 
the fact that \\[U, Pk]\\ = \\UP k U* — P k \\ allows one to estimate ||[C/, Pk]\\ above by 



max 



I <*j+i 1 2 M I Pj 1 2 - I Pj+i 1 2 1 } + max { I Pj<Xj - Pj+l<*j+l\}- 



18 



NATHANIAL P. BROWN 



But |a,f = |(1 - cos(^)), = |(1 + cos(^)) and = -| sin(^). Hence we 
get the inequality 

ll[tf,n]||<-^. 

Having the bounds on || [P fc , fi}\\ and ||[C/, PJH one completes the proof by a standard 
interpolation argument, similar to the proof of Theorem 7.1. □ 

Remark 8.3. We leave it to the reader to check that Pk(o~i) — > 5, h as k — > oo, for each ieZ. 
It will not always be the case that P k < Pfc+i but the fact that ||Pfc(t>) — f|| — * is good 
enough for all of the results in sections 3, 4 and 5 to hold (i.e. we really don't have to have 
an honest filtration in any of the proofs). 

Remark 8.4. Note that the rank of Pk is n k — m k and hence the numerical efficiency of this 
scheme will depend on how far out one must go in order to find long strings which are close 
together. 

A number of interesting examples naturally give rise to tri-diagonal operators on Z 2 (Z). 
For example, a natural discretization of ^ is -^{U — U*), where h is a fixed step size. A 
natural discretization of the one-dimensional differential operator (— A + c)/ = J^/ + cf, for 
some constant c, would be jjk(U + U*) + (c+ where / G Z°°(Z) is the constant function 
I(n) = 1 for all n. Also, Arveson's method of discretizing Schrodinger operators of one- 
dimensional quantum systems (cf. [1]) naturally leads to an analysis of operators of the form 
U + U* + / where f(n) = v(cos(n6)) for some continuous function v on the interval [—1,1] 
and some fixed positive number 6. Note that if - G Q then / will be a periodic function 
while in general it will be "almost" periodic and hence block balanced. Hopefully these 
examples justify spending the remainder of this paper analyzing how to implement Berg's 
technique on the class of tri-diagonal operators with periodic main diagonal and constant 
diagonals above and below the main diagonal. The interested reader should have no trouble 
now formulating a general result based on the previous lemma and the results of section 
4. However, while this method can be horribly inefficient in some cases, it works very well 
in the periodic tri-diagonal case. Moreover, for those interested in discretized Hamiltonians 
with ^ Q we will describe how to use continued fractions to find good "blocks" at the very 
end. (See also [5] for improved spectral computations in this case.) 

For the remainder of this paper T will be an operator of the form 

T = c\U + c^U* + /, 

where c±, c_i G C and / G Z°°(Z) is a periodic function with period p (i.e. f(n + p) — f(n), 
for all n G Z). Due to the periodicity of / and the constant coefficients on U and U* we may 
take Sk = 0, Sk = pk — 1, m k = — s fc and rik = 1 (notation as in the set-up for Lemma 8.2). 
Note that with these choices the "twisted" basis used in Berg's technique is 

Vfc = {ao :k d~m k + Po,kSi, a Sk ,kSo + (3 Sk:k o~ Sk }. 



QUASIDIAGONALITY 
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A tedious but straightforward computation shows that if P k is the orthogonal projection 
onto the span of V k then the matrix of P k TP k with respect to Berg's basis V k is 



/ o"o,o 


co,i 





■ ■ 





a 0,s k \ 


Cl,0 


^1,1 


01,2 


•• 











02,1 


02,2 


02,3 ■ • 








' 





03,2 


03,3 • • 

















■ ■ 


°s fc -l,s fc -l 












•• 


<T s k ,s k -l 


a s k ,s k J 



where 

a 3j = \®j,k\ 2 f( m k+j) + |/3j,fe| 2 /(l 

a j,j+± = { a j+i,k<Xj,k + Pj+i,k/3j,k)c-i, 
(8- 1 ) o-j+ij = (a^o^ + Pj,kPj+i,k)ci, 

O~0,s k = /9o,fc« Sfe ,fcCl, 
0s fe ,O = Po,kOlsk,kC-l- 

While it is unfortunate that the entries of the matrix P k TP k are more complicated than 
the entries of T it is a significant feature of Berg's technique that the property of being tri- 
diagonal (hence sparse) is almost preserved. Also note that \\[T, P k ]\\ < -|^- in this setting 
while the rank of P k is pk. Hence the error estimates are inversely proportional to the rank 
of the "almost" tri-diagonal matrix P k TP k . 

We end this paper by putting all the error estimates together in the case we are considering. 
The proof is immediate from Lemma 8.2 and the results of sections 3, 4 and 5. 

Theorem 8.5. Let T, P k e B{l 2 {Z)) be as above. 

(1) Let y G Z 2 (Z) be given and assume that T is invertible. If k > i(27r||T~ 1 || — 1) then 
the "almost" tri-diagonal matrix P k TP k is invertible and if x k is the unique solution 
to P k TP k x k = P k y then 

\\x-x k \\ < \\T- 1 \\(\\y-P k y\\ + ^^-), 

pk + 1 

where x is the unique solution to Tx = y. 

(2) For every e>0, a^(P k TP k ) -> a^{T) and, moreover, a {e \P k TP k ) C a {e+ ^(T). 

I iTy\\T\\ 

(3) a 2 (P k TP k ) - <7 2 (T) and a 2 (P k TP k ) cV wr a ^ T y 

2-7T 

(4) IfT — T* then a(P k TP k ) -> a(T) and a(P k TP k ) C^tt a (T). 

Remark 8.6. The theory of continued fractions gives a remarkable approximation property of 
irrational numbers (cf. [10]). Namely, if a G (0, 1) is an irrational number then there is a very 
simple and computable algorithm for finding coprime integers p n , q n such that la — — | < \. 

It follows that for any < j < q n we have the inequality \ja — j — \ < < — ■ Hence we 
can apply Berg's technique as above to a periodic approximation to a function of the form 
f{m) = cos(m9) and still have good control on the error estimates. 
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